Cyclotron Resonance in Strongly Magnetized Plasmas and Gamma Ray Burst 



S. Son 

18 Caleb Lane, Princeton, NJ 08540 



(N 
O 

> 
O 



C/3 . 



Sung Joon Moon 
28 Benjamin Rush Lane, Princeton, NJ 08540 
(Dated: November 6, 2012) 

A plausible scenario for the gamma ray and the hard x-ray burst in a strongly magnetized plasma, 
based on the collective plasma maser instability, is proposed. The physical parameters with which 
this scenario becomes relevant are estimated. The attractive feature of this scenario over the con- 
ventional cyclotron radiation theory is discussed. 
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INTRODUCTION 

According to a recent study, the thermal electron gyro- 
motion in a strongly magnetized plasma could lead to 
instabilities of electromagnetic (E&M) waves This 
instability, so-called the maser instability, has significant 
implications on strongly magnetized plasmas [2-LZ|. The 
goal of this paper is to estimate its relevance to the 
gamma ray or hard x-ray burst in the strongly magne- 
tized astrophysical plasmas [1, Q . The gamma ray burst 
in the absence of the magnetic field is not relevant to 
our analysis. Our analysis shows that the photons of 10 
keV to 1 MeV can be emitted from a relativistic plasma 
with the electron density 10^^ — 10^^ cm"'' and the mag- 
netic field of order of 10* to 10^'^ gauss. In contrast to the 
conventional cyclotron radiation where the photons are 
emitted rather uniformly, the angular distribution of the 
radiated photons could be sharply concentrated. Advan- 
tages of this possible phenomenon over the ones based on 
the incoherent cyclotron radiation are discussed. 

The rest of the paper is organized as follows. The in- 
stability growth rate is presented in Sec. 2, based on 
the Landau damping theory, which is analyzed when the 
photon wave vector is parallel to the magnetic field (Sec. 
3) and when the photon wave vector is not parallel to 
the magnetic field (Sec. 4). The gamma ray burst is dis- 
cussed in Sec. 5, and the difference between this theory 
and the conventional theories is presented Sec. 6, which 
is followed by summary (Sec. 7). 



COHERENT CYCLOTRON RADIATION 

Let us consider the growth rate of a collective E&M 
wave in the presence of a strong magnetic field Bqz; the 
detailed derivation when an E&M wave propagates paral- 
lel to the magnetic field is in Ref. [l[ . The physical mech- 
anism of the instability is similar to that of the maser or 
the free electron laser, where a coherent E&M wave gets 
amplified by the electron gyro-motion as manifested in 
the radio frequency wave of the solar corona 0, [l3| • The 



instability regime is ubiquitous in strongly magnetized 
plasmas ( [l| ) . To summarize briefly, the kinetic motion of 
a relativistic electron is solved in the presence of an E&M 
mode to the second order, and the ensemble average of 
the electron kinetic loss (gain) is taken over the electron 
distribution. Equating this energy loss to the E&M wave 
growth, the E&M growth (decay) rate is obtained. The 
growth rate for a general propagation direction 9, which 
is defined to be the angle between the magnetic field and 
the photon wave-vector, is 
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where k {lo) is the wave vector (frequency) of the E&M 
mode, fi = (/3i,/3||) = v/c, 0^ ^ pl + = 
ATTUf-e^ / rUe is the plasma frequency, S{f5) = cos(0) — 
Lo/ck -\- Wee /(cfc7 (/?)), ^ce = cBq / TTieC is the classical cy- 
clotron frequency, 7(/3) = (1 — (3'^)"'^^'^, f is the elec- 
tron distribution with the normalization of J fd^/S — 1, 
{A)s=o = J S{S)A(Pl3 is the integration in the velocity 
space with the constraint 5 = 0, and Qi and f22 are 
obtained below. Defining C to be the ratio of the wave 
energy density to the wave energy density in vacuum, ^ — 
Eiu/{E'^/Att), it is assumed in our analysis that C > 0, 
which is the case for most of E&M waves in dense plas- 
mas. For a generic angle 9, there could be two indepen- 
dent modes; TE and TM modes. The wave vector is given 
as k = k cos9z + sm9x. Let us define the TE (TM) mode 
as E = i?i cos(k ■ r — ujt)x and B = i?i {ck/u) cos(k ■ r — 
ujt) (cos 9y — sin 9z) (E = Ei cos(k • r — ojt) (cos 9x — sin 9z) 
and B = Ei{ck/uj) cos(k • r — iot)y). fli and il2 for the 
TE mode are 
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For the TM mode, they are given as 
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The non-relativistic limit of Eq. ([T]) is 
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where () is the ensemble average with — Vr ^ {itJ — 
u}ce)/k or (A) = J AS{v - Vr)cPv. When 6* = 0, the 
Eqs. ((T|and dU are reduced to the previously obtained 
results 

For more rigorous treatment when 6* 7^ 0, there are 
infinite expansions of the Bessel function J„(fc cos(0)rg), 
where Vg is the electron gyro-radius. However, the un- 
stable wave identified in our analysis has the condition 
kcos{9)rg < 1 so that Ji{k cos{9)rg) = 1. For this rea- 
son, Eq. ([T|) is derived assuming Ji = 1. 



CASE WHERE 9 = AND uj ^ ck 

Consider the case where the E&M field propagates par- 
allel with the magnetic field {9 = 0). There is no insta- 
bility if w = ck, because fli = 0. For simplicity, we 
refer to the first (second) term on the right-hand side of 
Eqs. HI) and ^ as the gyro-lasing (gyro-damping) term. 
For uj > cfc, the semi-classical instability growth rate is 
shown to be [H 
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where fithe = Vthe/c, Pr = Vr / c = {us - Wce)/c, Vthe IS the 
electron thermal velocity, and /|| {vr) — J f5{vz ~ Vr)(Pv. 
The instability criterion is F > 0. The maximum possible 
instability, ignoring the gyro-damping term, is estimated 
to be 
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where it is assumed that C — 1 [II- 

The Maxwellian distribution for fully relativistic elec- 
trons is f{f3) = exp(— 7A), where A = mpC^/T.. This is 
so-called the Maxwell- Juttner's distribution When 
A < 1, the distribution is peaked at 7 = 3/ A with 
the width of 6(3 = (A/3)^. We consider the case when 
to > ck > ujce- Assuming /3 = 1, the condition F > can 
be re-casted as 
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One necessary condition for the existence of the res- 
onance is given as w/cfc < Pmax + ^celln 



yj\ 4- uj'l^/c^k'^, where 7„ 



= (1 - Pla.)-"^ and 
If the electron distribu- 



= + (WceM)2 

tion has the relatively high slope near the resonance 
region {'^max — STe/nieC^), Eq. ([7]) is possible since 
l^pS] ^ 1 near the resonance. The growth rate without 
the gyro-damping term can be estimated, using V^/ = 
7^(3- A7)//? and /? ^ 1, as 
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where dS/dp = {P-VS)/\P\. 



CASE WHERE 6^0 AND uj = ck 

The resonance condition for the TM mode is = 
/3|| cos6'-t-a;ce/cfc7— 1 = 0. The gyro-lasing term vanishes 
at the resonance because fii = [c^k"^ /u)S cos{9) = 
from Eq. ([3]). However, the gyro-damping term does not 
vanish (f^a/cfc = cos6l(l - P\/2) - /3||). If < at 
5 = 0, this term acts as an amplifying term instead of 
a damping term. Consider the semi-classical case. From 
the resonance condition P^ cos 9 = l—uice/ck, Q.2 is given 
by 
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The instability condition [0,2 > 0) is given as P^ > 
(1 - Wee/ cfc) (1 - /3-L^/2). If P±=0, the condition is re- 

) .f\\{^r), (5) duced to Pr > \/l — LJce/ck, or COS0 = {l—UJce/ck)/Pr < 

/ a/1 —Li}^„ I ck. For a, ffiven Maxwellian distribution fnA 



\/^— uicejck. For a given Maxwellian distribution /m 
with the normalization J jdPz = 1 , the growth rate is 
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The growth rate has the maximum similar to Eq. (|6]) 
when Pthe < Pr < '^Pthe- For given Wee and cfc, the 
growth rate as a function of 9 has the maximum when 

Pthe COS = 1 - LOce/ck. 

For the TE mode, let us consider semi-classical elec- 
trons first. With the resonance condition Pr cos 9 = 
1 — ojee/cfc, we obtain Vli = ck {cos^ 9 — 1) and il2 — 
cfc(l - 131/2 ~ Pr cos 9) from Eq (g]). The instability cri- 
terion (F > 0) is 
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The maximum growth rate, ignoring the gyro-damping 
term, is given similar to Eq. ([5]). 

For the TE mode of the fully relativistic electrons, we 
consider only when 9 = 7r/2. The resonance surface is 
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given as S ■ 
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e/7w SO that uj = ujc 
Note also that fli 



,/7 and IV^S'I 
= —UJ and 



(1 - /3j/2)cfc from Eq. ©. A shmlar analysis to Eq. © 
can be used, and the instability criterion (F > 1) becomes 
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Note that V^gS-V^/ < for a positive gyro-lasing term 
since fii < 0. Assuming (3 = 1, Eq. ([T^ and is simplified 
to 
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In the case of the Maxwellian plasma, \^pfM\/fM = 
(37^ - A73)/3, where A = m^c^ /T^. If A7 < 2, Eq. ^ 
and \7pS ■ Vjsf < are satisfied. Assuming /3 = 1, it is 
estimated that 
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Note that df /d(3 can be as large as 7"^, and the maximum 
growth rate is given as 
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The comparison between the case when 9 = and 
when 7^ is in order. For the case when 9 = 0, ft2 
vanishes ii uj = ck, thus the instability can occur only 
when UJ ^ ck; this requires a very high electron density 
for gamma rays or hard x-rays. On the other hands, 
when 9 ^ 0, an explosive instability is possible even when 
UJ = ck ^ ujpe- Therefore, the case when 9 ^ Q is more 
probable than the case when 9 = 0. 



GAMMA RAY BURST 

The maximum growth rate of the instability in the 
classical plasma is given in Eq. ([6]), and the one in the 
fully relativistic plasmas is in Eqs. dH) and ([15]). For non- 
Maxwellian plasmas, the growth rate should be derived 
based on Eqs. ([T|) and (|3]). In this section, it is assumed 
that there is sufficient free energy in the plasma so that 
the concern on the Gardner's constraint is not relevant [1, 

Sill- 
Regardless of whether the plasma is Maxwellian or not, 
the maximum instability growth rate can be summarized 
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where 5(7) = 1. The energy loss rate via the cyclotron 
radiation is given as 
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where P is the loss power and Td is the ratio of the loss 
rate to the electron perpendicular kinetic energy. The 
condition for the instability growth rate to exceed that 
of the cyclotron radiation is 



(18) 



where Tmax is given in Eq. ^ or (|T6|). Note that the 
left-hand side of Eq. (|18|) is proportional to the electron 
density, while the other side is not; the instability dom- 
inates the cyclotron radiation for higher densities. For 
the non-relativistic electrons, the ratio Vmax/^d is given 
as = 1.79 X (fie/k^). For the relativistic electrons, it is 
given as Tmax/Td = iieg{l)/k^l'^- 

Let us define the reference frame Sq where the electron 
average drift is zero. Let us assume that the resonant 
velocity Vr — {c — ujce)/k <iC c in the 5o-frame. For a 
non-relativistic Maxwellian plasma in the S'o-frame, the 
density satisfying the ratio Tmax/^d = 1 for 0.01 keV 
(10 keV) photon is = 10^^ cm"^ (io25 cm"3) r^^^ 
photon with the energy hck in the S'o-frame is observed as 
a photon with the energy ^ohck in the observer's frame, 
where the 70 is the Lorentz factor between the S'o-frame 
and the observer's frame. Noting that the relativistic 
factor 7o is usually between 10 and 1000 according to 
the general literature on the gamma ray burst the 
photon of energy 10 keV < johck < 1 MeV in the 
observer's frame corresponds to the photons of energy 
0.01 keV < 7ofi.cA; < 100 keV in the So-frame, whose 
cyclotron frequency corresponds to the magnetic field of 
10^ to 10" gauss. 



THE INSTABILITY THEORY VERSUS THE 
INCOHERENT CYCLOTRON RADIATION 
THEORY 

The theoretical consideration in our study has the fol- 
lowing distinctive features, compared to the conventional 
cyclotron radiation theory: 

{i) The fastest growing E&M mode would be the one 
parallel (or perpendicular) to the magnetic field. The 
intensity of a collective E&M mode is proportional to 
exp{jtcos9) when 9 = and to exp(7icos(6' — 7r/2)) 
when 9 = 7r/2, so that the peak angle is narrow when 
7t ^ 1. Therefore, in a certain direction, more intense 
photons would be observed than predicted from the con- 
ventional incoherent cyclotron radiation, whose intensity 
is proportional to cos^ 9. The actual energy power re- 
quirement for the gamma ray burst might be less than 
that of the conventional theory Q . 

(m) As the electrons emit the photons via the instabil- 
ity, the temperature becomes further anisotropic, render- 
ing the E&M wave perpendicular to the magnetic field 
unstable to the Weibel instability. In turn, the Weibel in- 
stability mitigates the temperature anisotropy, emitting 
low-frequency photons in the perpendicular direction. 



4 



[iii) Our analysis suggests that the gamma ray burst 
could be originated from a more compact and dense ob- 
ject than conventionally believed. In a highly dense 
plasma considered here, an incoherent photon has a short 
mean free path due to the gyro-damping as well as the 
Thomson scattering, while the collective photons could 
overcome these dampings through the gyro-lasing. In 
other words, the compact and dense plasma could be op- 
tically thin for the coherent photons, while optically thick 
for the incoherent photons. 

To be more specific, let us consider the Thomson scat- 
tering. The decay rate is given as Tt — UpaTC = 
ne^gC, where = /mc.c^ is the classical electron ra- 
dius. From Eq. (fT6l) . the ratio Fmax/rT is given as 



Tmax/rT — i^/kre). Then, the Thomson scattering is 
weaker than the gyro-lasing if the photon wavelength is 
longer than the classical electron radius, whose photon 
energy roughly corresponds to 100 MeV, which is the case 
in our scenario. This suggests that the photon can get 
amplified, overcoming the Thomson scattering when the 
gyro-lasing from the instability is strong. However, as the 
instability gets weaker, the Thomson scattering becomes 
eventually the dominant decay mechanism of the pho- 
tons, where the photon mean free path is / = (l/ngr^). 
For example, for Uf, = 10^^/cc, I = 10^ cm. For a given 
object, the collective photons could be excited inside the 
object and sustained via the gyro-lasing. But, if there is 
no instability on the surface of the object, the photons 
will decay due to the Thomson scattering within rather 
short mean free path; in order for the photons to escape 
from the object, it is necessary that the gyro-lasing from 
the instability exists on the surface of the object. The 
delicate interplay between the Thomson scattering and 
the gyro-lasing in the process of the collective photons 
escaping from an object is beyond the scope of this pa- 
per. 

{iv) The growth rate of the instability is sensitive to 
the slope of the distribution function at the resonance. 
Consider a shock region where two plasmas of different 
drifts violently encounter. It is plausible that the parallel 
and perpendicular electron temperatures are comparable 
but the distribution could have two sharp humps in the 
parallel direction, at which case the instability could be 
strong. 

(u) Finally, our theory provides an escape scenario of 
photons from dense plasmas in a changing magnetic field. 
A typical example would be the decreasing magnetic field 
(e.g., the foot point of the solar corona). Consider the 
TE mode when a photon propagates with 9 = tt /2 from 
r = and the magnetic field decreasing with the increas- 
ing r (or B = B{r)z with B{ri) < B{r2) if ri > r2)- 
From the resonance condition u — ujceHi the relativistic 
factor 7 of a resonant electron should decrease with the 
decreasing Wce- Assuming the electron temperature re- 
mains the same along the photon propagation, the gyro- 
lasing term stays positive if A7 < 2 at r = (as analyzed 



in Sec. 4). Then, collective photons are radiated primar- 
ily in the direction of (perpendicular to) the magnetic 
field as photons from the incoherent cyclotron radiation 
suffer a severe damping and photons not parallel with (or 
perpendicular to) the magnetic field experience a weaker 
gyro-lasing term than the photons parallel (or perpendic- 
ular) to the magnetic field. 



SUMMARY 

A scenario of the gamma ray burst, based on the recent 
radiation theory JL] , is proposed and examined. The pre- 
vious analysis [l| is generalized to an arbitrary angle. The 
estimation shows that the coherent burst of 10 keV to 1 
MeV photons is plausible in a relativistic plasma when 
7o = 10 ~ 1000 and the electron density is higher than 
some critical value, Ue > lO'^® ~ 10^® cm~^. It is shown 
that a rather compact dense object with less available 
energy could cause the short gamma ray burst and that 
the observed gamma rays would be coherent rather than 
incoherent. In addition, the coherent photons of lower 
frequency comparable to the plasma frequency might be 
observed due to the Weibel instability. This is partic- 
ularly relevant to the case when 6 — ti/2 because the 
low frequency coherent photons (high frequency coherent 
photons) from the Weibel instability (from the instability 
studied here) can be observed simultaneously. The above 
conjecture might be useful in verifying whether the sce- 
nario proposed here would account for some of the short 
gamma ray burst events observed in the satellites 

While our estimation is rather focused on the gamma 
and hard x-rays, a similar mechanism would manifest in 
generating soft x-rays in the inertial confinement fusion 
plasma IJ]. The electron beam of 7 > 10 ~ 100 and the 
magnetic field of 10^ gauss can be readily generated in 
laboratories. Even a magnetic field of 10^ gauss might 
be possible [lH. Then, the photon generated from the 
instability may have energy between 10 eV and 1 keV. 
Complications would be t he elect ron quantum diffraction 
effect and the degeneracy 
is in progress. 
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